Methods and systems for location estimation

ABSTRACT

Methods and systems for location estimation. Signals are received from a plurality of satellites and a plurality of base stations by a mobile station. The location of the mobile station is estimated according to TOA (Time-Of-Arrival) and TDOA (Time Difference-Of-Arrival) measurements calculated based on the signals.

BACKGROUND

The present disclosure relates generally to methods and systems for location estimation, and more particularly, to methods and systems that estimate locations of mobile stations based on network-based and satellite-based signals.

Wireless location technologies have drawn a significant amount of attention over the past few decades. Different types of Location-Based Services (LBSs) have been proposed, including emergency 911 (E-911) subscriber safety services, location-based billing, navigation system, and applications for the Intelligent Transportation System (ITS). However, to meet all applications' needs, it is required to provide more accurate location estimation of a mobile device, such as a mobile phone, a smart phone, and a PDA under different environments.

A variety of wireless location techniques have been proposed. These techniques can be classified into network-based and satellite-based location schemes. The representative algorithms for the network-based location techniques are TDOA (Time Difference-Of-Arrival) and AOA (Angle-Of-Arrival) schemes. FIG. 1 is a schematic diagram illustrating a network-based architecture. In the TDOA scheme, a mobile station (MS) 110 receives signals coming from different wireless base stations (121, 122, and 123), and measures the time difference of the signals. In the AOA scheme, a home base station of the mobile station 110 observes the arrival angle of signals from the mobile station 110. On the other hand, the well-adapted technology for the satellite-based location method is GPS (Global Positioning System). FIG. 2 is a schematic diagram illustrating a satellite-based architecture. The mobile station 110 receives signals coming from a plurality of satellites (131 and 132), and measures TOA (Time-Of-Arrival) of the signals.

The performance of the location techniques, however, varies under different environments. For example, at rural areas, less base stations in a network-based system are installed, resulting in weak incoming signals or shortage of signal sources. Therefore, at rural areas, the network-based (i.e. TDOA, AOA) methods result in degraded performance for the location determination of mobile devices. On the other hand, the major problem for the satellite-based (i.e. GPS) method is that the performance considerably degrades when the satellite signals are severely blocked in urban valley areas.

SUMMARY

Therefore, it is desired to have a better location estimation method and system that adapt all kinds of environments. In the present invention, methods and systems for hybrid location estimation of mobile devices are provided.

In an exemplary embodiment of a method for location estimation, signals from a plurality of satellites are received by a mobile station, and TOA measurements are performed accordingly. The first position of the mobile station is estimated using an expansion method based on the TOA measurements. Signals from a plurality of base stations are received by the mobile station, and TDOA measurements are performed accordingly. The second position of the mobile station is estimated using the expansion method based on the TDOA measurements. The first and second positions are incorporated to estimate the location of the mobile station.

In an exemplary embodiment of a method for location estimation, signals from a plurality of satellites and a plurality of base stations are received by a mobile station. The location of the mobile station is estimated accordingly.

An exemplary embodiment of a system for location estimation comprises a plurality of satellites, a plurality of base stations, and a mobile station. The mobile station receives signals from the satellites, performs TOA measurements accordingly, and estimates the first position of the mobile station using an expansion method based on the TOA measurements. The mobile station further receives signals from the base stations, performs TDOA measurements accordingly, and estimates the second position of the mobile station using the expansion method based on the TDOA measurements. The mobile station incorporates the first and second positions to estimate the location thereof.

An exemplary embodiment of a system for location estimation comprises a plurality of satellites, a plurality of base stations, and a mobile station. The mobile station receives signals from the satellites, performs TOA measurements accordingly, and estimates a first position of the mobile station using an expansion method based on the TOA measurements. The mobile station further receives signals from the base stations, performs TDOA measurements accordingly, and transmits the first position of the mobile station and the TDOA measurements to a specific base station. The specific base station estimates a second position of the mobile station using the expansion method based on the TDOA measurements, and incorporating the first and second positions to estimate the location of the mobile station.

Methods for location estimation may take the form of program code embodied in tangible media. When the program code is loaded into and executed by a machine, the machine becomes an apparatus for practicing the disclosed method.

DESCRIPTION OF THE DRAWINGS

The invention will become more fully understood by referring to the following detailed description with reference to the accompanying drawings, wherein:

FIG. 1 is a schematic diagram illustrating a network-based architecture;

FIG. 2 is a schematic diagram illustrating a satellite-based architecture;

FIG. 3 is a schematic diagram illustrating an embodiment of a system for location estimation;

FIG. 4 is a flowchart of an embodiment of a method for location estimation;

FIG. 5 is a schematic diagram illustrating an embodiment of a MS-assisted system for location estimation;

FIG. 6 is a schematic diagram illustrating an embodiment of a MS-based system for location estimation; and

FIG. 7 is a schematic diagram illustrating the relationship between a mobile station and a home base station.

DESCRIPTION

Methods and systems for location estimation are provided.

FIG. 3 is a schematic diagram illustrating an embodiment of a system for location estimation. The system combines the signals from both the GPS system and the cellular networks.

The system comprises a mobile station (MS) 310, such as a mobile phone, a smart phone, and a PDA, a plurality of base stations (BSs) (321, 322, and 323), and a plurality of satellites (331 and 332). The MS 310 is equipped with a GPS receiver (not shown). Additionally, each BS has a forward-link pilot channel that continuously broadcasts its pilot signals to provide timing and phase information for the MS 310 in the cellular network. Each BS has a dedicated reverse-link pilot channel from the MS 310 to provide initial acquisition, time tracking, and power control measurements. Each BS is equipped with antenna arrays for adaptive beam steering to facilitate the AOA measurement. The TDOA measurements are conducted at the MS 310 by obtaining the signals via the forward-link pilot channels from the BSs. It is called mobile-based mode. It should be noted that, depending on the capability of mobile devices, the TDOA measurements can also be conducted at the BSs, naming the mobile-assisted mode. On the other hand, the AOA signals are transmitted from the MS 310 to the BSs using the respective reverse-link pilot channels. The AOA measurements are performed at respective BS using its antenna arrays for two-dimensional adaptive beam steering. In some embodiments, to avoid signal degradation due to the near-far effect, the home BS of the MS 310 provides the capability of the AOA measurement, but is not limited thereto.

FIG. 4 is a flowchart of an embodiment of a method for location estimation. The GPS-equipped MS receives signals from the satellites and conducts TOA measurements (t_(k)) (401). The MS processes the TOA measurements and estimates the three-dimensional position (i.e. X_(G)=[x_(G) y_(G) z_(G)]^(T)) of the MS using an expansion method (402). In this embodiment, the two-step least square (LS) method is employed, but is not limited thereto. Other approaches, such as using Taylor Expansion, Range-Scale Algorithm, Believable Factor Algorithm, and LLOP, may be adopted for the estimation. Additionally, TDOA measurements (t_(i,j)) are performed at the MS by obtaining signals from its home BS and the neighboring BSs via the forward-link pilot channels (403). The AOA measurement (θ) is conducted at the home BS by receiving the signals from the MS via the reverse-link channel (404). The TDOA and AOA measurements are combined to estimate three-dimensional position (X_(C)) of the MS using the two-step LS method (405). A filtering technique is performed on the respective estimated positions (X_(G) and X_(C)), to smooth measurement noise and track position data from both the TOA and the TDOA/AOA channels (406 and 407). In this embodiment, a Kalman filtering technique is employed, but is not limited thereto. Data fusion is then performed to incorporate both the means ({overscore ({circumflex over (X)})}_(G)) and ({overscore ({circumflex over (X)})}_(C)) of the filtered estimations ({circumflex over (X)}_(G) and {circumflex over (X)}_(C)) from the TOA and the TDOA/AOA measurements based on their signal variations (408). The fused position estimate ({circumflex over (X)}_(f)) of the MS is obtained (409).

The location estimation of MS can be performed in either MS-assisted or MS-based mode. The two modes of systems depend on the requirement of the communication bandwidth and the computation power of the MS. The location estimation of the present invention can be applied to either MS-assisted or MS-based systems.

FIG. 5 is a schematic diagram illustrating an embodiment of a MS-assisted system for location estimation. This type of architecture is suitable for MS with insufficient computation capability. The GPS-equipped MS receives signals from the satellites and conducts TOA measurements (t_(k)). The MS processes the TOA measurements and estimates the three-dimensional position (X_(G)) of the MS using a two-step LS method. Additionally, the TDOA measurements are calculated at the MS by obtaining signals from its home BS and the neighboring BSs via the forward-link pilot channels. These two sets of information, the location estimation (X_(G)) from the GPS system and the TDOA measurements (t_(i,j)) from the cellular network, are transmitted back to the home BS via the reverse-link pilot channel. The AOA measurement (θ) is conducted at the home BS by receiving the signals from the MS via the reverse-link channel. A location server (not shown) at the home BS performs location estimation by combining the AOA and the TDOA measurements. The two-step LS method is utilized to estimate three-dimensional position (X_(C)) of the MS. The location server performs a Kalman filtering technique to smooth measurement noise and track position data from both the TOA and the TDOA/AOA channels, and performs data fusion to incorporate both the means ({overscore ({circumflex over (X)})}_(G)) and ({overscore ({circumflex over (X)})}_(C)) of the filtered estimations ({circumflex over (X)}_(G) and {circumflex over (X)}_(C)) from the TOA and the TDOA/AOA measurements based on their signal variations. The fused position estimate ({circumflex over (X)}_(f)) of the MS is obtained.

FIG. 6 is a schematic diagram illustrating an embodiment of a MS-based system for location estimation. This type of architecture is suitable for MS that possesses adequate computation capability. The AOA measurement (θ) is obtained from the home BS and is transmitted to the MS via the forward link pilot channel. From the forward-link pilot channel, the MS computes the TDOA measurements (t_(i,j)) based on signals from its home and neighboring BSs. The MS performs location estimate (X_(C)) by combing the TDOA and the AOA measurements using the two-step LS method. The GPS-equipped MS receives signals from the satellites and provides location estimate (X_(G)) based on the TOA measurements (t_(k)). The Kalman filtering techniques are performed for both the TOA and the TDOA/AOA channels at the MS to obtain location estimates, {circumflex over (X)}_(G) and {circumflex over (X)}_(C). The final position estimation ({circumflex over (X)}_(f)) is acquired after the fusion results from means of the location estimates, {overscore ({circumflex over (X)})}_(G) and {overscore ({circumflex over (X)})}_(C).

The details of TOA, TDOA and AOA measurements, location estimation using the two-step LS method, Kalman filtering, and data fusion follow.

TOA, TDOA, and AOA Measurements

The mathematical models for the TOA, TDOA, and AOA measurements are presented. The three dimensional (3-D) coordinates of the MS are exploited in the proposed hybrid location estimation algorithm.

The 3-D TOA measurement t_(k) from the GPS system is obtained by $\begin{matrix} {{t_{k} = {{\frac{1}{c}r_{k,o}} + n_{k}}},{k = 1},2,\ldots\quad,N,} & (1) \end{matrix}$ where c is the speed of light, the subscript k represents the measurement with regard to the kth satellite, the subscript o represents the original time, and n_(k) is the measurement noise associated with the TOA measurement t_(k). The relative distance (r_(k,o)) between the MS and the kth satellite can be obtained as r _(k,o)=√{square root over ((x−x _(k))²+(y−y _(k))²+(z−z _(k))²)}  (2), where (x,y,z) represents the MS position, and (x_(k),y_(k),z_(k)) is the location of the kth satellite.

The 3-D cellular-based TODA measurement t_(i,j) is obtained by computing the time difference between the MS with regard to the ith and the jth BSs: $\begin{matrix} {{t_{i,j} = {{\frac{1}{c}\left( {r_{i,o} - r_{j,o}} \right)} + n_{i,j}}},} & (3) \end{matrix}$ where i and j represent the ith and the jth BS, and n_(i,j) is the measurement noise of t_(i,j).

In some embodiments, the antenna arrays at the home BS can measure the x and y directions, the AOA measurement (θ) of the cellular system is represented as $\begin{matrix} {{\theta = {{\tan^{- 1}\left( \frac{y - y_{1}}{x - x_{1}} \right)} + n_{\theta}}},} & (4) \end{matrix}$ where θ is the horizontal angle between the MS and its home BS, (x₁,y₁) is the horizontal coordinate of the home BS, and n_(θ) is the measurement noise of θ. Location Estimation Using Two-Step LS Method 3-D TOA Location Estimation

In some embodiments, it is required to obtain signals from at least four satellites (N≧4) to solve for the two-step LS problem from the TOA measurement. The 3-D TOA measurement as in (1) can be rewritten as Hx=J  (5), where ${H = \begin{bmatrix} {{- 2}x_{1}} & {{- 2}y_{1}} & {{- 2}z_{1}} & 1 \\ {{- 2}x_{2}} & {{- 2}y_{2}} & {{- 2}z_{2}} & 1 \\ \vdots & \vdots & \vdots & \vdots \\ {{- 2}x_{N}} & {{- 2}y_{N}} & {{- 2}z_{N}} & 1 \end{bmatrix}},{x = \begin{bmatrix} x & y & z & R \end{bmatrix}^{T}},{J = \begin{bmatrix} {r_{1}^{2} - L_{1}} \\ {r_{2}^{2} - L_{2}} \\ \vdots \\ {r_{N}^{2} - L_{N}} \end{bmatrix}},$ R=x²+y²+z², L_(k)=x_(k) ²+y_(k) ²+z_(k) ², and r_(k)=ct_(k) for k=1,2, . . . ,N. Where r_(k) is the measured distance between the MS and the kth satellite.

The concept of the two-step LS method is to acquire an intermediate location estimate in the first step, assuming that x, y, z and R are not correlated. The second step of the method releases this assumption by adjusting the intermediate result to obtain an improved location estimate, X_(G)=[x_(G) y_(G) z_(G)]^(T). Further details of the two-step LS method are omitted here.

3-D TDOA/AOA Location Estimation

In order to evaluate the feasibility of using the 2-D AOA measurement for location determination of MS, it is necessary to investigate the geometric relationship between the MS and the home BS. FIG. 7 is a schematic diagram illustrating the relationship between a mobile station and a home base station. As shown in FIG. 7, the position of the home BS is defined as the reference point, and the geometric relationship resulting from AOA measurement noise n_(θ) _(2D) in x-y plane is shown as: D ₁ sin n _(θ) _(2D) =(x−x ₁)sin θ−(y−y ₁)cos θ  (6), where D₁=√{square root over ((x−x₁)²+(y−y₁)²)}. r_(1,o) represents the relative distance between the MS and the home BS. When D₁>>h, D₁≈r_(1,o). Hence, the approximation is D₁ sin n_(θ) _(2D) =r_(1,o) sin n_(θ) _(2D)   (7).

Therefore, (6) is changed as r _(1,o) sin n _(θ) _(2D) ≈(x−x ₁)sin θ−(y−y ₁)cos θ  (8).

Due to sin n_(θ) _(2D) ≈n_(θ) _(2D) when |n_(θ) _(2D) |<<1, (8) can be approximately rewritten in linear form as 0≈−(x−x ₁)sin θ+(y−y ₁)cos θ+r _(1,o) n _(θ) _(2D)   (9). The benefit of using the AOA measurement from the home BS is to reduce the requirement to obtain signals from one extra BS. Due to the limitation of the AOA measurement obtained from the two-dimensional antenna array, the vertical angle between the MS and the home BS, however, is not measurable. The undetermined vertical angle may result in excess errors in the geometric relationship (9), especially where the MS is close to the home BS location. Therefore, two cases are considered. In some embodiments, the AOA measurement further contains a distance-determining device (not shown), which determines whether the long distance case or the short distance case should be applied. a) Long Distance between the Home BS and the MS

In this embodiment, at least three TDOA measurements from the BSs and one AOA measurement from the home BS are required for location estimation. The 3-D TDOA and AOA measurements from (3) and (4) can be reformulated as (5), where ${H = {- \begin{bmatrix} {x_{2} - x_{1}} & {y_{2} - y_{1}} & {z_{2} - z_{1}} & r_{2,1} \\ {x_{3} - x_{1}} & {y_{3} - y_{1}} & {z_{3} - z_{1}} & r_{3,1} \\ \vdots & \vdots & \vdots & \vdots \\ {x_{N} - x_{1}} & {y_{N} - y_{1}} & {z_{N} - z_{1}} & r_{N,1} \\ {{- \sin}\quad\theta} & {\cos\quad\theta} & 0 & 0 \end{bmatrix}}},{x = \begin{bmatrix} x & y & z & r_{1,o} \end{bmatrix}^{T}},{and}$ ${J = {\frac{1}{2}\begin{bmatrix} {r_{2,1}^{2} - L_{2} + L_{1}} \\ {r_{3,1}^{2} - L_{3} + L_{1}} \\ \vdots \\ {r_{N,1}^{2} - L_{N} + L_{1}} \\ {{2x_{1}\sin\quad\theta} - {2y_{1}\cos\quad\theta}} \end{bmatrix}}},$ and r_(i,1)=ct_(i,1) for i=2,3, . . . , N. r_(1,o) represents the relative distance between the MS and the home BS as in (2). b) Short Distance between the Home BS and the MS

Since it is not feasible to utilize the AOA measurement when the MS is in the neighborhood of the home BS, an additional 3-D TDOA measurement is required for location estimation. Hence, in this embodiment, at least four TDOA measurements are required for location estimation. The 3-D TDOA measurements form (2) can be written as (5) where ${H = {- \begin{bmatrix} {x_{2} - x_{1}} & {y_{2} - y_{1}} & {z_{2} - z_{1}} & r_{2,1} \\ {x_{3} - x_{1}} & {y_{3} - y_{1}} & {z_{3} - z_{1}} & r_{3,1} \\ \vdots & \vdots & \vdots & \vdots \\ {x_{N} - x_{1}} & {y_{N} - y_{1}} & {z_{N} - z_{1}} & r_{N,1} \end{bmatrix}}},{x = \begin{bmatrix} x & y & z & r_{1,o} \end{bmatrix}^{T}},{and}$ $J = {{\frac{1}{2}\begin{bmatrix} {r_{2,1}^{2} - L_{2} + L_{1}} \\ {r_{3,1}^{2} - L_{3} + L_{1}} \\ \vdots \\ {r_{N,1}^{2} - L_{N} + L_{1}} \end{bmatrix}}.}$ Kalman Filtering

The Kalman filtering technique is employed in this embodiment for post-processing of signals. It provides the capabilities of range measurement, smoothing, and noise mitigation for the TOA and TDOA/AOA data. Furthermore, the Kalman filter can also be utilized for trajectory tracking of the MS. The measurement process and the state equation for the Kalman filter can be written as X _(k) =C{circumflex over (X)} _(k) +e _(k)  (10), {circumflex over (X)} _(k) =A{circumflex over (X)} _(k-1) +w _(k-1)  (11), where X_(k)=[x_(k) y_(k) z_(k)]^(T) is the measured data from either {circumflex over (X)}_(C) or {circumflex over (X)}_(G) at the time instant k, {circumflex over (X)}_(k)=[{circumflex over (x)}_(k) ŷ_(k) {circumflex over (z)}_(k) v_(x) _(k) v_(y) _(k) v_(z) _(k) ]^(T) is the estimated state vector, and e_(k) and w_(k) are the process and measurement noises. The matrix C in (10) relates the state ({circumflex over (X)}_(k)) to the measurement (X_(k)) as $C = {\begin{bmatrix} I & {\Delta\quad{tI}} \\ 0 & I \end{bmatrix}.}$ The matrix A in (11) is the state transition matrix as M=[I 0]. This embodiment applies the Kalman filtering technique to both TOA and TDOA/AOA channels, but is not limited thereto. Data Fusion

The main function of data fusion is to merge disparate types of information to enhance positioning accuracy. By merging the TOA and TDOA/AOA estimates, the resulting data provides feasible location estimates under different environments (e.g. urban, suburban, and rural). In this embodiment, the fusion process is based on the Bayesian Inference model, which improves an estimate with known signal variations. It is understood that the fusion process is not limited to the Bayesian Inference model. The final position estimate and its variance are obtained as follows: $\begin{matrix} {{\hat{X}}_{f} = {{\frac{\sigma_{G}^{2}}{\sigma_{G}^{2} + \sigma_{C}^{2}}{\hat{\overset{\_}{X}}}_{C}} + {\frac{\sigma_{C}^{2}}{\sigma_{G}^{2} + \sigma_{C}^{2}}{\hat{\overset{\_}{X}}}_{G}}}} & (12) \\ {and} & \quad \\ {{\sigma_{f}^{2} = \frac{1}{\frac{1}{\sigma_{C}^{2}} + \frac{1}{\sigma_{G}^{2}}}},} & (13) \end{matrix}$ where {overscore ({circumflex over (X)})}_(C)=[{overscore ({circumflex over (x)})}_(C) {overscore (ŷ_(C) {overscore ({circumflex over (z)})}_(C)]^(T) is the mean of the location estimate {circumflex over (X)}_(C) from the TDOA/AOA channels, {overscore ({circumflex over (X)})}_(G)=[{overscore ({circumflex over (x)})}_(G) {overscore (ŷ_(G) {overscore ({circumflex over (z)})}_(G)]^(T) is the mean of the location estimate {circumflex over (X)}_(G) for the TOA channel, σ_(C) is the standard deviation of {circumflex over (X)}_(C), and σ_(G) is the standard deviation of {circumflex over (X)}_(G).

The hybrid location estimation determines the MS location by combining outcomes from both network-based and satellite-based techniques. In addition to the longitude and latitude of MS, the altitude information of MS is also obtainable therefrom. The accuracy of location estimation can thus be improved.

Methods for determining operational settings, or certain aspects or portions thereof, may take the form of program code (i.e., executable instructions) embodied in tangible media, such as products, floppy diskettes, CD-ROMS, hard drives, or any other machine-readable storage medium, wherein, when the program code is loaded into and executed by a machine, such as a computer, the machine thereby becomes an apparatus for practicing the methods. The methods may also be embodied in the form of program code transmitted over some transmission medium, such as electrical wiring or cabling, through fiber optics, or via any other form of transmission, wherein, when the program code is received and loaded into and executed by a machine, such as a computer, the machine becomes an apparatus for practicing the disclosed methods. When implemented on a general-purpose processor, the program code combines with the processor to provide a unique apparatus that operates analogously to application specific logic circuits.

While the invention has been described by way of example and in terms of preferred embodiment, it is to be understood that the invention is not limited thereto. Those skilled in the technology can still make various alterations and modifications without departing from the scope and spirit of this invention. Therefore, the scope of the present invention shall be defined and protected by the following claims and their equivalents. 

1. A method for location estimation of a mobile device, comprising: performing TOA (Time-Of-Arrival) measurements according to signals from a plurality of satellites; estimating a first position using an expansion method based on the TOA measurements; performing TDOA (Time Difference-Of-Arrival) measurements according to signals from a plurality of base stations; estimating a second position using the expansion method based on the TDOA measurements; and incorporating the first and second positions to estimate the location of the mobile station.
 2. The method of claim 1 further comprising: smoothing the first and second positions by a filtering technique.
 3. The method of claim 2 wherein the filtering technique is Kalman filtering.
 4. The method of claim 1 further comprising: performing an AOA (Angle-Of-Arrival) measurement by a specific base station; and estimating the second position using the expansion method based on the TDOA and AOA measurements.
 5. The method of claim 4 wherein the specific base station is a home base station of the mobile device.
 6. The method of claim 4 further comprising: determining whether the distance between the mobile device and the specific base station is above a threshold.
 7. The method of claim 1 further comprising: estimating the first position using the expansion method based on the TOA measurements and the locations of respective satellites.
 8. The method of claim 1 further comprising: estimating the second position using the expansion method based on the TDOA measurements and the locations of respective base stations.
 9. The method of claim 1 wherein the expansion method is two-step Least Square (LS) method.
 10. The method of claim 1 wherein the first and second positions are incorporated by a Bayestian statistical inference model to estimate the location of the mobile station.
 11. A method for location estimation of a mobile station, comprising: receiving signals from a plurality of satellites and a plurality of base stations; and estimating the location of the mobile station.
 12. The method of claim 11 further comprising: performing TOA (Time-Of-Arrival) measurements according to the signals from a plurality of satellites; performing TDOA (Time Difference-Of-Arrival) measurements according to the signals from a plurality of base stations; and estimating the location of the mobile station according to the TOA and TDOA measurements.
 13. The method of claim 12 wherein the TOA and TDOA measurements are incorporated by a Bayestian statistical inference model to estimate the location of the mobile station.
 14. The method of claim 12 further comprising: performing an AOA (Angle-Of-Arrival) measurement by a specific base station; and estimating the location of the mobile station according to the TOA, TDOA, and AOA measurements.
 15. The method of claim 14 wherein the TOA, TDOA, and AOA measurements are incorporated by a Bayestian statistical inference model to estimate the location of the mobile station.
 16. The method of claim 14 further comprising: determining whether the distance between the mobile station and the specific base station is above a threshold.
 17. The method of claim 14 wherein the specific base station is a home base station of the mobile device.
 18. The method of claim 12 further comprising: smoothing the estimation according to the TOA and TDOA measurements by a filtering technique.
 19. The method of claim 18 wherein the filtering technique is Kalman filtering.
 20. The method of claim 14 further comprising: smoothing the estimation according to the TOA, TDOA and AOA measurements by a filtering technique.
 21. The method of claim 20 wherein the filtering technique is Kalman filtering.
 22. The method of claim 12 further comprising: estimating the location of the mobile station using an expansion method based on the TOA and TDOA measurements and the locations of respective satellites and base stations.
 23. The method of claim 22 wherein the expansion method is two-step Least Square (LS) method.
 24. The method of claim 14 further comprising: estimating the location of the mobile station using an expansion method based on the TOA, TDOA, and AOA measurements and the locations of respective satellites and base stations.
 25. The method of claim 24 wherein the expansion method is two-step Least Square (LS) method.
 26. A system for location estimation, comprising: a plurality of satellites; a plurality of base stations; and a mobile station, performing TOA (Time-Of-Arrival) measurements according to signals from the satellites, estimating a first position using an expansion method based on the TOA measurements, performing TDOA (Time Difference-Of-Arrival) measurements according to signals from the base stations, estimating a second position using the expansion method based on the TDOA measurements, and incorporating the first and second positions to estimate the location thereof.
 27. The system of claim 26 wherein the mobile station further performs a filtering technique on the first and second positions, respectively, thus to smooth measurement noise therefrom.
 28. The system of claim 27 wherein the filtering technique is Kalman filtering.
 29. The system of claim 26 wherein a specific base station further receives signals sent from the mobile station, performs an AOA (Angle-Of-Arrival) measurement accordingly, and transmits the AOA measurement to the mobile station.
 30. The system of claim 29 wherein the specific base station is a home base station of the mobile station.
 31. The system of claim 29 wherein one of the base stations further determines whether the distance between the mobile device and the specific base station is above a threshold.
 32. The system of claim 29 wherein the mobile station further estimates the second position using the expansion method based on the TDOA and AOA measurements.
 33. The system of claim 26 wherein the mobile station estimates the first position using the expansion method based on the TOA measurements and the locations of respective satellites.
 34. The system of claim 26 wherein the mobile station estimates the second position using the expansion method based on the TDOA measurements and the locations of respective base stations.
 35. The system of claim 26 wherein the expansion method is two-step Least Square (LS) method.
 36. The system of claim 26 wherein the first and second positions are incorporated by a Bayestian statistical inference model to estimate the location of the mobile station.
 37. A system for location estimation, comprising: a plurality of satellites; a plurality of base stations; and a mobile station, performing TOA (Time-Of-Arrival) measurements according to signals from the satellites, estimating a first position using a two-step Least Square (LS) method based on the TOA measurements, performing TDOA (Time Difference-Of-Arrival) measurements according to signals from the base stations, and transmitting the first position and the TDOA measurements to a specific base station, wherein the specific base station estimates a second position using the two-step Least Square (LS) method based on the TDOA measurements, and incorporates the first and second positions to estimate the location of the mobile station.
 38. The system of claim 37 wherein the specific base station is a home base station of the mobile station.
 39. The system of claim 37 wherein the specific base station further smoothes the first and second positions by a Kalman filtering method.
 40. The system of claim 37 wherein the specific base station further, performs an AOA (Angle-Of-Arrival) measurement according to signals from the mobile station, and estimates the second position using the two-step Least Square (LS) method based on the TDOA and AOA measurements.
 41. The system of claim 37 wherein the specific base station incorporates the first and second positions by a Bayestian statistical inference model.
 42. A device, comprising: means for performing TOA (Time-Of-Arrival) measurements according to signals from a plurality of satellites; means for estimating a first position based on the TOA measurements; means for performing TDOA (Time Difference-Of-Arrival) measurements according to signals from a plurality of base stations; means for estimating a second position based on the TDOA measurements; and means for incorporating the first and second positions to estimate the location of the device.
 43. The device of claim 42 wherein the incorporating means incorporates the first and second positions by a Bayestian statistical inference model.
 44. A device, comprising: means for performing TOA (Time-Of-Arrival) measurements according to signals from a plurality of satellites; means for estimating a first position based on the TOA measurements; means for performing TDOA (Time Difference-Of-Arrival) and AOA (Angle-Of-Arrival) measurements according to signals from a plurality of base stations; and wherein a home base station of the device estimates a second position based on the TDOA and AOA measurements, and incorporates the first and second positions to estimate the location of the device.
 45. A device, comprising means for estimating the location of the device according to signals from a plurality of satellites and a plurality of base stations. 